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1 The intrinsic complexity of chemical models 

Photodissociation region (PDR) models are used to understand the evolution of the 
FUV illuminated matter both in our Galaxy and in external galaxies. To prepare for 
the unprecedented spatial and spectroscopic capabilities of ALMA and Herschel, two 
different kinds of progresses are currently taking place in the field. First, numerical 
models describing the chemistry of a molecular cloud are being benchmarked against 
each other to ensure that all models agree not only qualitatively but also quantitatively 
on at least simple case fQ (RoUig et al. 2007). Second, new or improved chemical 
rates are being calculated/measured by several theoretical and experimental groups. 
However, the difficulty of this last effort implies that only a few reactions (among the 
thousand ones used in chemical networks) can be thoroughly studied. New numerical 
tools are thus being developed for taking into account the impact of the uncertainties 
of the chemical rates on the chemical model predictions, and their comparison with 
observed abundances [e.g. Wakelam et al. 2005, 2006). In view of the intrinsic 
complexity of building reliable chemical networks and models, there is an obvious need 
of well-defined observations that can serve as basic references. PDRs are particularly 
well suited to serve as references because they make the link between diffuse and dark 
clouds, thus enabling to probe a large variety of physical and chemical processes. 

2 The Horsehead edge as a chemical laboratory 

The illuminated edge (PDR) of the western condensation presents one of the sharpest 
infrared filament (width: 10" or 0.02 pc) detected in our Galaxy by ISOCAM. The 
most straightforward explanation given by Abergel et al. (2003) is that most of the 
dense material is within a flat structure viewed edge-on and illuminated in the plane of 
the sky by cOri. The H2 fluorescent emission observed by Habart et al. (2005) is even 
sharper (width: 5"), implying the inclination of the PDR on the plane-of-sky to be 
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Fig. 1. Emission maps obtained with the IRAM Plateau de Bure Interferometer or 30m 
single-dish, except for the H2 v=l-0 S(l) emission observed with the NTT/SOFI and the 
PAH mid-IR emission observed with IS0-LW2. The maps have been rotated by 14° counter- 
clockwise around the image center to bring the exciting star direction in the horizontal 
direction as this eases the comparison of the PDR tracer stratifications. Maps have also 
been horizontally shifted by 20" (compared to the projection center: RA — 05h40m54.27s, 
Dec = -02°28'00", J2000) to set the horizontal zero at the PDR edge delineated as the 
vertical red line. Either the synthesized beam or the single dish beam is plotted in the 
bottom left corner. The emission of all the lines is integrated between 10.1 and 11.1 kms^^. 
Values of contour level are shown on each image lookup table (contours of the H2 image 
have been computed on an image smoothed to 5" resolution). The green crosses display the 
positions where we derived column densities and abundances (see Table [l]). 



less than 5°. The Horsehead ridge thus offers the opportunity to study at small linear 
scales (1" corresponds to 0.002 pc at 400 pc) the physics and chemistry of a PDR with 
a simple geometry, very close to the prototypical kind of source needed to serve as a 
reference to chemical models. Since 2001, we started to study the Horsehead PDR 
mainly with the IRAM Plateau de Bure interferometer at 3 mm and the IRAM-30m at 
1mm achieving spatial resolutions from 3 to 11" (except the H2 fluorescent emission 
observed at l"-resolution with the NTT/SOFI instrument). Figure [T] displays all 
the observed, high resolution maps already acquired. Those maps trace the different 
layered structures predicted by photochemical models according to chemical reaction 
networks, excitation conditions and radiative transfer. 

Abergel et al. (2003) deduced from the distance between crOri and the PDR that 
the intensity of the incident far UV radiation field is x ~ 60 relative to the interstellar 
radiation field in Drainc's units. Through the modelling of the H2 and CO emission, 
Habart et al. (2005) showed that the PDR has a very steep density gradient, rising 
to riH ^ 10^ cm~'^ in less than 10" {i.e. 0.02 pc), at a roughly constant pressure of 
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Table 1. Column densities and abundances of several chemical species from the UV- 
illuminated PDR (upper part) to the shielded region (lower part) of the Horsehead edge. 
Coordinates offsets are given both in the Equatorial system from RA = 05h40m54.27s, Dec 
= -02°28'00" (J2000) and in the coordinate system adapted to the source geometry and used 
in Figure [1] Abundances are computed with respect to the density of protons (nn). 



P ~ 4 X 10^ Kcm"'^. These observations were followed by a chemical study of small 
hydrocarbons (CCH, C-C3H2, C4H). Pety et al. (2005) showed that the abundances 
of the hydrocarbons are higher than the predictions based on pure gas phase chemical 
models (Le Petit et al. 2002 and reference therein). These results could be explained 
either by the photoerosion of the large aromatic molecules and/or small carbon grains 
(C. Joblin, 2007, this volume) or by a turbulent mixing which would transport in the 
illuminated part of the PDR molecules after their productions in the dark part (M. 
Gerin, 2007, this volume). Goicoechea et al. (2006) linked the PDR model predictions 
with detailed non-LTE, nonlocal excitation and radiative transfer models adapted 
to the Horsehead geometry. They showed that the gas sulfur depletion invoked to 
account for CS and HCS"*" abundances is orders of magnitude lower than in previous 
studies of the sulfur chemistry. 

Finally, Pety et al. (2007) studied the deuterium fractionation in the Horsehead 
edge from observations of several H^^CO+ and DCO+ lines. A large [DCO+]/[HCO+] 
abundance ratio (> 0.02) is inferred at the DCO+ emission peak, a condensation 
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shielded from the illuminating far-UV radiation field where the gas must be cold 
(10-20 K) and dense (nn > 4 x lO^cm"^). DCO"*" is not detected in the warmer 
photodissociation front, implying a lower [DCO''"]/[HCO+] ratio (< 10^^). To our 
knowledge, this is the brightest DCO"'" emission (4 K) detected in an interstellar cloud 
close (angular distance < 40") to a bright H2/PAH emitting region. This opens the 
interesting possibility to probe at high resolution the chemical transition from far-UV 
photodominated gas to "dark cloud" shielded gas in a small field of view. 

3 Towards an observational benchmark 

An ideal observational benchmark would deliver to chemists a set of abundances 
(with the associated uncertainties) as a function of the distance (or extinction) to 
the illuminating star. This goal is difficult to achieve for several reasons: 1) The 
geometry of the source is never as simple as wished when it is known at all; 2) 
The spectra produced by the instruments must be inverted to obtain abundances; 3) 
The spectra are often measured at very different angular resolutions, implying beam 
dilution and/or mixing of different gas components. For several years now, we have 
started to systematically study the western edge of the Horsehead nebula because its 
geometry is not only well understood but also quite simple (almost ID and viewed 
edge-on). The density profile across the PDR is well constrained and there are several 
current efforts to constrain the thermal profile. The combination of low distance to 
Earth (400 pc), low illumination (x 60) and high density {n ~ lO^cm"'^) implies 
that all the interesting physical and chemical processes can be probed in a field-of-view 
of less than 50" with typical spatial scales ranging between 1 and 10". 

All those observations are done by the same team, using the same instruments 
(mainly IRAM Plateau de Bure and 30m) and the same methods both of data reduc- 
tion and data analysis. For each species, we are trying to observe several transitions 
at similar angular resolutions (from 5 to 15") to constrain properly the excitation 
conditions and derive accurate column densities and abundances. Obtaining emission 
maps when possible has proved essential to understand the spatial distributions of 
the species. We are preparing the public release of cuts of spectra so that the radia- 
tive transfer analysis can be refined as knowledge of coUisional rates progresses. In 
the meantime. Table [1] is our first attempt at quantitatively summarizing the results 
obtained up to now. In the future, the zoom capacity and resolving power of ALMA 
will enable to measure many different specie transitions at a resolution of 1", enabling 
to resolve all the physical and chemical gradients. 
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